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configurations the role of dislocation multipoles cannot be neglected, but the theory 
presented can still be applied. 
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1. Introduction 

The macroscopic plastic deformation of a crystalline material is the result of the 
movement of a huge number of dislocations interacting with each other through the 
long-range anisotropic shear stress field generated by a single dislocation. This feature 
and the fact that dislocation movement is often constrained to single glide plane result 
in dynamics exhibiting very high spatiotemporal complexity. This is indicated, for 
instance, by the multitude of observed dislocation patterns or by the recently revealed 
fact that plastic deformation is characterised by intermittent strain bursts with scale-free 
size distribution [1, 2, 3, 4]. 

Statistical descriptions of dislocation systems mainly concentrate on deriving 
continuum models for the evolution of the different dislocation densities. Several such 
theories have been proposed so far both for two [5, 6, 7, 8, 9] and three dimensions 
[10, 11, 12, 13, 14]. In all these models a key problem is to take into account 
dislocation-dislocation correlations during the coarse graining of the microstructure. 
(It is known that neglecting the dislocation-dislocation correlations leads to unphysical 
results [8,9].) It is found that correlation effects enter into the theories through gradient- 
like stress terms [6, 8, 14], which were already proposed earlier phcnomenologically 
[15, 16]. Motivated by its crucial role many investigations on the properties of dislocation 
correlations have been published recently [17, 18, 19, 20, 21]. 

Beside the one-particle dislocation density function the probability distribution 
function of internal stresses serves as another important quantity for the statistical 
characterisation of dislocation systems. Previously this function was determined for 
relaxed 2D dislocation configurations [22, 23] as well as for fractal-hke dislocation 
morphologies [24], but in these studies the role of external shear stress was not 
investigated. For many aspects, however, the external stress obviously plays a major 
role. 

What gives the stress distribution function outstanding importance is the fact that 
it can be directly measured on real materials by X-ray diffraction, even in situ. Briefly, 
it was shown by Groma and co-workers that the shape of the broadened Bragg peaks 
depend on the distribution of the internal strain [25, 26, 27] that is proportional to 
the stress. Consequently, from the Bragg peaks the stress distribution function can be 
extracted. 

Another practical issue is the numerical modelling of discrete dislocations. It was 
observed by Groma and Bako that the force acting on a dislocation can be divided 
into a slowly varying component coming from the further dislocations and a stochastic 
contribution of the near dislocations [28]. In the proposed 0{N) stochastic dislocation 
dynamics simulation method the latter is drawn from the stress distribution function 
[22, 29]. 

Like in the previous studies on the topic [22, 23, 24] we consider two dimensional 
(2D) dislocation systems. This makes the calculations much simpler, but according to 
X-ray studies, under quite general conditions the results are applicable for deformed 
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3D crystals as well. On the other hand, there is a growing interest towards recently 
observed physical systems with real 2D underlying lattices whose macroscopic properties 
are highly influenced by dislocations. Examples include dusty plasma crystals where 
particles are arranged in a hexagonal lattice [30, 31], vortex lattices in type II 
superconducting films [32], colloidal crystals [33, 34] and foams [35]. 

In this paper the stress distribution function is investigated in the presence of 
external loading. In section 2, theoretical calculations are carried out to arrive at 
a closed analytical expression for the Fourier transform of the distribution function. 
This expression depends on the spatial dislocation-dislocation correlation function. 
Analytical expression has not been found yet for the correlation function, it can be 
determined only numerically [17, 21]. To arrive at close form for the Fourier transform 
of the stress distribution function, like in an earlier study [23], in section 3 we consider 
randomly distributed monodisperse dislocation dipole systems. This permits to give 
the asymptotic decay of the stress distribution (section 4). It is found that due to 
external load the stress distribution function becomes asymmetric. Namely, a l/(r|r|^) 
hke term proportional with the apphed external stress is added to the 1/|t|^ hke tail of 
the distribution function (r denotes the shear stress). In the vicinity of the origin the 
distribution function is shifted with a value which is again proportional to the applied 
external stress. 

Stress distribution functions obtained numerically are presented in section 5. 
First, the theoretical predictions are validated on the monodisperse dipole systems. 
Then investigations of relaxed configurations obtained by discrete dislocation dynamics 
simulations follow. It is found that the asymptotic decay seen on monodisperse systems 
remains valid. Difference is found only in the numerical values of some parameters, 
which is attributed to the large number of dislocation multipoles present in the relaxed 
configurations making the system more resistant against the external stress. 

2. The analytical form of the stress distribution function 

The general mathematical formulation of the stress distribution function was given 
previously by Groma and Bako [22]. In this section, first, we briefiy present their 
results then continue with analysis specific to the problem subject of this paper. 

2.1. General description 

Let us consider a system of parallel straight edge dislocations with line directions 
parallel to the z axis. For the sake of simphcity we restrict our analysis for single glide, 
so the Burgers vector of the dislocations can be only ±6 with b taken to be parallel 
to the X axis. Under these conditions, the dislocations' movement is parallel with the 
X axis and the three dimensional problem can be treated as a 2D one in the xy plane. 
The position and the sign of the Burgers vector of the ith dislocation in the xy plane is 
denoted by r^, and Si (r^ e M^, Sj e {—1, 1} and i e {1, . . . , N}), respectively. 
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Assuming linear elasticity the internal shear stress field generated by the N 
dislocations at a given point r e is 

TV 

<V-''^(r, ri, . . . , Vn) := ^ Sin^^ir - r^), (1) 

i=l 

where Tind is the shear stress field in the xy plane generated by an individual positive 
sign edge dislocation positioned in the origin: 

^ . ~ ^^ sin(y^)cos(2y^) 

rind(r-) := Gb ^^2 + ^2)2 = ' ' (2) 

in which G is a combination of the elastic moduli: G — , where /i and v are the 

shear modulus and Poisson's ratio, respectively. It has to be noted, that the given form 
of Tind is valid only for dislocations in an infinite medium. 

The subject of this paper is to determine the stress distribution function Pstr at 
a randomly chosen point r. By definition Pstr(T, r)dr is the probability of finding the 
stress generated by the system between the limits 

T - Y < <V-''"(r-, n, . . . , rjv) < T + y. (3) 



To determine Pgtr Markoff 's method [36] is followed. 

Let w^j^''"'^'^ denote the N particle distribution function of the system. After 
introducing the function 

A. - dx) :=( - ^ <<' ■■""('■•'■'•■■■•'■«)< ^ + ^- . 

I 0, otherwise 

(4) 

Pstr can be expressed in the form 

Pstr(T,r)dT= Yl y"dVi...y"dV^A^i'-'^^(r,ri,...,r^,T,dT) 

xw'^'-''-{r^,...,r^). (5) 

It was previously pointed out by Groma and Bako that, in contrast to Markoff's 
original deduction, the dislocation-dislocation correlations must be taken into account 
to avoid the system size dependence of the stress distribution function [22]. Then the 
resulting form of the Fourier transform of the probability density function {P^j.[q, r) 
^/M^str(r,r)e-^^dT) is 

ln(Ps;(?,r))= - 5] /■pf(n)S-(r-n,g)dVi + ^ 5] J] J j p\^{r,)p\\r,) 
X ^^'^^(ri, r2)S^i(r - n, q)B'^{r - r^, g)dVidV2 + . . . , (6) 

where 

P'(r, g) := 1 - expOgsTindlr)), (7) 
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and dV^'^'^ denotes the dislocation-dislocation correlation functions defined as 



' ^ ^' pr(n)pr(r,) ' 

in which and P2^'*^ are the one and the two particle dislocation density functions, 
respectively. The superscripts in these functions refer to the sign of the dislocations, 
e.g., {r-i, is proportional to the probability of finding a positive dislocation at ri 
and a negative one at r2- 

In order to evaluate (6), further assumptions have to be made. In the rest of 
this paper we consider only neutral infinite homogeneous configurations, meaning p\ is 
constant in space: 

P%.s---P{{r) (9) 
and the densities of the positive and negative sign dislocations are equal: 

Pdis := 2p+, = 2p-,,. (10) 

Consequently, 

• The direct r dependence of Pgtr must vanish and 

• The (^2^'*^ correlation functions depend only on the difference of their arguments 
[17, 21]: 

dr\r,,r2)^dr'{r,-r2). (11) 

We note that for symmetry reasons d2~{r) = d2~^{r) and d2~^{r) = d2~{—r) hold 
for all r values, thus the correlations can be described using only the rf^^ and the 
d2~ functions [17, 21]. As in previous studies, we neglect the three and higher order 
correlations meaning that the terms not written out explicitly in (6) are omitted. 

2.2. The symmetric part of the distribution's Fourier transform 

In the absence of external stress the dislocation pair correlation functions are symmetric: 
d2'''\r) = dl^''\-r) with si,S2 e {+,-} [17, 21]. Then, according to (6), P^f, is real, 
and therefore the distribution function Pgtr is symmetric. This case was thoroughly 
studied recently by Csikor and Groma [23]. Starting from (6), for relaxed dislocation 
configurations they found asymptotes for P^^ in the limits of j^j — > and \q\ oo. 

Under applied external shear stress the correlation functions of the system become 
different (for details see section 3.1) which modifies the second term in (6). After a 
short calculation one concludes that the real part of (6) depends only on and 
the combination c?^' + d2^ . Since these functions do not change considerably due to 
the external stress (see section 3.1 for details), we assume that Re (Pj;) can be well 
approximated by the form derived by Csikor and Groma 



Re[ln(P,f,(5))] = < 



Cpdisg^ In ( — ) , if |g| ^ 0, 
-yPdisl^l, if|g|^oo. 



(12) 
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where C — j{Gb)^. D is a parameter which can only be determined numerically (for 
relaxed systems D = 1.35G6pJ°'^ was obtained) and geflf is a constant [23] but it does 
not play role in the further considerations. 



2.3. The antisymmetric part of the distribution's Fourier transform 



After substituting defined by (7) into (6) and performing a few straightforward 
transformations one obtains 



Im [\n{P^l(q))]^(p^^/2f j d+-(d)T(d,q)d% 



where 



T{d,q) 



:=/sin(, 



d 



Tind r - 



- Tind r + 



d 



dV, 



(13) 



(14) 



in which it is assumed that (i^"^(r) = d^'^(— r), and so the second term in (6) vanishes 
for si = S2 = 1 and Si = S2 = — 1- This assumption will be verified in section 3.1. 

The main subject of this paper is the evaluation of (13). If no external stress 
is applied, the right hand side of (13) vanishes [23], but in a loaded system this 
term plays an important role. In the following, first, the behaviour of function T is 
investigated by performing the integral in (14) numerically. Afterwards, the properties 
of Im [in (-Pstr(?))] analysed in detail. 

Since analytical solution for the integral appearing in (14) was not found, it is 
calculated numerically. Because of the unusual form of the integral, however, the 
numerical integration is not trivial. First, we have to make an important remark. 

Integration in (14) should be performed on the whole plane, but with any 
numerical algorithm one can calculate the integral for only a finite region and study its 
region size dependence. Formally speaking, the integral of an arbitrary / : ^ M 
function is approximated as: 

R 

/(r)d^r = ^im J drr J dipf{r,ip). (15) 

R2 [0,2n[ 

This approximation inherently assumes that the / function can be integrated 
successively (first for if then for r) which is not true for all possible integrablc / functions. 
According to Fubini's theorem, for successive integration, among others. 



/ 



r|/(r, (p)\dr < +oo 



(16) 



must hold [37]. 

Since r^d exhibits a type decay, the argument of the sine function in (14) has 
a type asymptote as r — > oo. This implies that for the first term of the expansion 
of the sine function in (14), which also has a tail, the condition given by (16) is 
not fulfilled. So, according to the aforementioned theorem, the integral in (14) cannot 
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be performed successively, and therefore it cannot be determined by a straightforward 
numerical algorithm: 



T{d,q) 



sm I q 



d 



d 



dV 



27r 



^ind 



d 

r 

2 



d 

^+2 



(17) 



7^ J drr J d</7sin 



In order to overcome this problem we have to note that 

• The higher order terms in the expansion of the sine function in (14) have r" 
. . . type singularities and so can be integrated successively. 

• The integral of the first term [which breaks the condition (16)] must be zero, because 
Tind is an odd function, and hence 



-6 ^-10 



d 



Tind r - 



- Tind r + 



d 



dV 



M2 



d^r — q Ti 



' ind 



r + ^ ) d^r 



0. 



Accordingly, if the first term of the series of the sine function is subtracted from 
the argument of the integral in (14), then what remains can already be integrated 
successively. So, 



T{d,q) 



sm q 



d 



d 



Tind I r- + - I - Tind I r - - 



dV 



OO 2-7T 




sm q 



Tind 



- Tind r - 



d 







-Q 



Tind 



d\ f d 



r dip dr. 



(19) 



2 J V 2 

With this, the determination of the function T is now simplified to the numerical 
integration of the right hand side of (19). During the numerical computation, however, 
special attention is needed at the points r = ±| since the argument of the sine function 
has a 1/r type singularities there. Moreover, the argument exhibits a strong angular 
dependence around ±|, too. As a result, the first term of the integrand in (19) oscillates 
rapidly between 1 and —1 near the mentioned points. To handle this behaviour correctly, 
close to the points r = ±| the integration grid-point distance was adaptively reduced. 
Formally, denoting by Ax the distance of the neighbouring grid-points at a given grid- 
point if the condition 



q/S.x ■ V, 



Tind 



d 



r — 



< 



2tt 



(20) 



is not fulfilled additional grid-points are introduced. [Equation (20) corresponds to the 
criterion that there must be at least A^b grid-points in every period of the oscillating 
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function.] The refinement of the mesh has to be stopped at a certain distance from 
the singular points. However, it is easy to see that the integral for a small symmetric 
domain around ±| vanishes as the size of that domain approaches zero. The mesh used 
is demonstrated in figure 1. 

1 

0.5 

^ 
-0.5 

-1 

-1 -0.5 0.5 1 

X 

Figure 1. The grid-points of the algorithm developed for evaluating the integral in 

/ Q 55 \ 

(19) next to the points d = [ ' ] . Closer grid-points are taken near the singular 

V 0-5 J 

points r = ±d/2. 




The function T obtained by the method explained above is plotted in figure 2. For 
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40 60 
q<Gb) 



80 100 



Figure 2. The T function obtained numerically at d 
\q\ ~* oo a linear curve was fitted. 



0.4 
0.5 



In the limit of 



the two asymptotic regimes T can be well described by the forms: 



T{d,q) 



{Ghf ■a{d)q^\Yi 



Gh ■ (3{d)q, 



if |g| ^ 0, 
if |g| oo. 



(21) 
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Table 1. The values of the a, qo and /? functions at different d-s obtained by fitting 
to numeric data. 







a{d) 


On (d) 


did) 


0.4 


0.5 


-0.334 


0.181 


0.819 


0.425 


0.5 


-0.248 


0.132 


0.861 


0.45 


0.5 


-0.f63 


0.0685 


0.902 


0.475 


0.5 


-0.0807 


8.90 • 10-3 


0.942 


0.5 


0.5 


a 


a 


0.981 


0.525 


0.5 


0.075 


46.5 


1.021 


0.55 


0.5 


0.f47 


5.55 


1.059 


0.575 


0.5 


0.214 


2.88 


1.097 


0.6 


0.5 


0.276 


2.12 


1.135 



^ The values are not given due to the huge numerical errors. 

According to figure 2 the linear relation in (21) at |g| — ^ oo needs no explanation, 
while the Iq'I one is verified in figure 3, where T{d,q)/q^ is plotted as a function 
of ln(g). The obtained straight line proves (21). By fitting the above functions to the 
numerical data of T the values of the parameters a{d), /3(d) and qo{d) can be determined 
for an arbitrary d. In table 1 values for these three parameters at different d values 
corresponding to dipole angles close to 45° are displayed. (The reason for investigating 
only the vicinity of 45° will be discussed in detail in section 3.) 

q-Gb 

0.01 0.1 1 




-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 
Iniq-Gb) 



Figure 3. T{d,q)/ (Gbq^) is plotted as a function of h\(q) at d = ( ' ). The 

V 0-5 / 

resulting linear function confirms the \q\ — > asymptote in (21). 

The following scaling relation can be easily derived directly from the definition of 
the T function given by (14): 

T{k-d,q) = k'-T(^d,^y (22) 
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From this one can find immediately, that according to (21) the relations 

a{k ■ d) = — 

qo{k ■d) = k- qo{d), 
(3{k-d) =k-(3{d) 

hold for an arbitrary A; G M. Furthermore, as a consequence of (14) 

T{d,q) = -T{d,q) 



(23) 
(24) 
(25) 

(26) 



is fulfilled too, where d 



-4 



According to the relations (23)-(26) in order to determine the complete form of 
functions a, (3 and go, one has to investigate their values only at dipole angles between 
and 90°. Using the data from table 1, the following approximate functions were found 
(see figure 4): 



-a' 




Figure 4. The numerically obtained values of the functions a, go and /3, and the fitted 
functions given by (27)- (29). 



a((ia:,0.5) 



141 



1.70 [-exp(-1.80(|4| - 0.5)) + 1.00] 



(27) 
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/3(4,0.5) 



_4 

14 



(1.58- 141 + 0.191). 



(29) 



From the relations (23) and (27) one gets 

a(4, dy) = a (^^4, 2rfj^ ■ 0.5^ = (^-^, 0.5^ 



■ \dy 



■0.85 



- exp -0.90 



dv 



- 1.0 



1.00 



The functions Qq and 6 can be obtained on a similar way leading to 

/ , , N , , , , /0.218 ■ \dy\\ 
qo[dx, dy) = 1.23 • \dy\ ■ exp (^"jj-pq^ J , 



P{dx,dy) 



dx 
14 



;i.58- 141 + 0.38- \dy\). 



(30) 

(31) 

(32) 



By combining (21), (30), (31) and (32) one concludes 

{Gb) 



Im(ln(PX(5))) 



3 (Text) 3, / Qoi'^cxt] 

PdisQ in 



Go PdisQ: 



if |g| 



if |g| — > oo. 



(33) 



where 



«'(re.t) := ^ / 4-(d)a(d)dd, 



ext j 



/ dt~{d)a{d) \n{qo{d))dd 

R2 

J dt'{d)a{d)dd 



and 



/3'(Text) 



Pdis 

2 



4-(d)/3(d)d(i. 



(34) 



(35) 



(36) 



M2 



3. The effect of the external stress on the microstructure 

To arrive at the form of the Fourier transform of the stress distribution function that is 
comparable with the ones obtained on real relaxed dislocation systems the values of the 
parameters a', q'^ and (3' have to be determined. According to their definitions given 
by (34)-(36), they depend on the correlation function So, in order to get a closed 
form for the probability distribution P^^., the properties of ci^~ and the microstructure 
of relaxed dislocation systems have to be analysed. 
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3.1. The pair correlation function 

The (i^^ correlation function of a relaxed homogeneous dislocation system can be 
directly determined by discrete dislocation dynamics simulations by counting the 
relative coordinates between positive and negative dislocations in the equilibrium 
configuration obtained numerically. In order to get smooth correlation function one 
has to perform averaging over many different realizations. Figure 5(a) shows the 
correlation function at zero external stress. In figure 5(b) d^" at external shear stress 
Text = (l/\/128) ■ Gby/pdis ~ 0.09 ■ Gby/pdis is plotted. In both cases the simulations 
were started from a random distribution of 64 positive and 64 negative dislocations and 
2000 different realizations were used for averaging. 




10 



0.1 



-3-2-10123 

X [Pdis ] 



(b) 



1 






1 







10 



0.1 



-3-2-10123 

X [Pdis J 



1 + (4 + 4"')/2 




-3-2-10123 

X [Pdis ] 



10 



0.1 



Figure 5. The 1 + correlation functions obtained numerically from many 
different relaxed dislocation configurations, (a) at zero external stress (b) at Text = 
(1/^128) • Gh^/p^s ~ 0.09 • Gh^/p^s external shear stress, (c) 1 + (dj" + d^+)/2 
is also plotted at Text = (l/-\/128) • Gh^p^is ~ 0.09 • Gh^fp^ to demonstrate, that 
it practically does not change due to the external stress. In all cases averaging was 
performed over 2000 different realizations. 



One can see in figure 5 that in relaxed systems the positive and negative dislocations 
tend to form narrow dislocation dipoles. Applied external stress makes asymmetric 
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but its nature close to the origin remains unchanged. The reason for this is that sphtting 
a narrow dipole requires very high stress. Even when the acting stress is a bit below 
the yield stress, only wide dipoles are split up (in the 2D single slip dislocation systems 
under study the yield stress is Xy 0.1 • Gb^/p^ [38]). We note that the symmetric part 
of [which is equal to (c?^~ + ^2 ^)/2] is nearly unaffected by the external stress (see 
figure 5(c)). 

The d^"*" correlation function was also plotted in figure 6 to demonstrate that it 
practically does not change due to external stress, and that dislocation walls are also 
present in stressed configurations. In figure 7 a stressed and an unstressed relaxed 
configuration is shown. The presence of the mentioned dislocation dipoles and walls can 
be clearly seen. 



(a) 




X [Pdis ] 



10 



0.1 



(b) 




-1 1 

^ [Pdis ] 



10 



0.1 



Figure 6. The 1 + rff correlation functions obtained numerically from 2000 different 
relaxed dislocation configurations. Plot (a) and (b) correspond to the same stress level 
as the two in figure 5. 



3.2. Monodisperse dislocation dipole systems 

The properties of the ci^^ correlation function mentioned above and the relaxed 
configurations plotted in figure 7 indicates that mainly dislocation dipoles affect the 
values of a', and both in stressed and unstressed systems [see equations (34)- 
(36)]. So, in this section as a first step we consider a set of randomly positioned 
ideal 45° dipoles with same momenta. This will be referred to as a monodisperse 
dislocation dipole system. It is important to note, that this is only a model system, 
since it is not in a dynamic equilibrium. As a result of random positions, the dipoles 
may even overlap. Although, this approach may seem to be oversimplified at first 
glance, as it is demonstrated later the results obtained can capture the properties of 
the stress distribution function found on more realistic configurations. Similar method 
was used earlier by Csikor and Groma [23]. In their paper the relaxed 2D distribution 
was envisaged as a set of ideal dislocation dipoles and short ideal walls. Using this 



(b) 
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Figure 7. Relaxed dislocation configurations obtained by discrete dislocation 
simulations. In figure (a) there was no applied shear stress, while in figure (b) the 
external shear stress was Text = (1/V128) • Gby/pdis ~ 0.09 • Gby/pdis- Notice that in 
the stressed system the narrow dipoles are still present as well as the dislocation walls. 



concept the form of the distribution function was given analytically for the stress free 
case. The resulting theoretical distribution function was in complete agreement with 
the one obtained numerically. Since small external stress does not split up the dipoles, 
it is reasonable to adopt this approach for the loaded case considered in this paper. 

We start our analysis by studying the behaviour of an isolated dipole subjected 
to shear stress. In equilibrium the sum of the forces (the Peach-Koehler force and the 
force produced by the external stress) acting on both dislocations must be zero. If 
there is no external stress there are four possible stable dislocation dipole orientations 
corresponding to the dipole angles 45°, 135°, 225° and 315°. If external stress is applied 
let e^(Text,l/o) denote the change of the x coordinate of the negative dislocation (see 
figure 8) with yo denoting the distance of the shp planes of the dislocations. 

The e^(Tcxt,yo) functions can be determined from (2) (see figure 9). It is easy to 
see that 

s^(rext,l/o)^-2^ if|rext|-0. (37) 

It should be noted that the domain of e^{., yo) is the interval 

stresses the dipole sphts. 

The monodisperse dipole system under external stress is built up of randomly 

positioned, separately deformed dipoles, as in figure 9. The main advantage of 

investigating a monodisperse system is that the corresponding c?^~ correlation function 

can be exactly given. If there are N/2 positive and N/2 negative dislocations then 

_ 2 2 
dt'ir) = 6{y - yo)6{x - y - e+{T^^t,yo)) - tt, (38) 

Pdis 

where ?/o > 0. 



Gh Gb 
iyo ' 4j/o 



since at higher 
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Figure 8. Under external shear stress Text the stress free dipole moment d becomes 
d'. The change of the x coordinate is denoted by e~(Tcxt, yo) for (a) 135° and (c) 225° 
dipoles and by s~^{Tc^t,yo) for (b) 45° and (d) 315° dipoles, where yo = dy. 



-0.5 













e (-JoVyo 






-2Text>'0/(Gfo) 













-0.3 -0.2 -0.1 0.1 0.2 0.3 
^extVCG^') 

Figure 9. The functions. Their tangent in the origin given by (37) is also plotted. 



From (30), (38) and (34) we arrive at 



X 



fI ■ \y\ 



rO.85 



0.85 / 6{y~yo)- 



M2 

exp (^-0.90 

— exp 



X 

y 

0.90 

\y\ 



1.0 -h 1.00 



da; dy 



£+(rcxt,yo) +1.00 



dy 
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0.85 

— exp 

\yo\ 

0.85-0.90 



-0.90 

\yo\ 



£+(rext,l/o) +1.00 



|yo| 



£ (Text,yo) ~ -1-53 



(39) 



/3'(Text) = ^J 4"(r-)/?(r)dr = J S{y - yo)5{x - y - £+(Te.t,yo)) 



X (1.58- \x\ 

\x\ 



0.38- 



da; dy 



= 1.58 • [yo + £+(Text, yo)] + 0.38 • yo ~ 1-96 • yo - 3.16^yo'. (40) 

The above results correspond to configurations consisting of only 45° dipoles. After a 
short and straightforward calculation one finds that for systems consisting of the four 
possible ideal dipoles in equal number: 

.Text 



Q;'(Text) ~ -1.53 



Gh 



P'ir^^O ^ -3.16^yo^ 



(41) 
(42) 



(The function could be determined in a similar way but it will not play any role in 
the further investigations so, its actual form is not given explicitly.) 



3.3. Disperse dislocation dipole systems 

A more realistic model of a 2D relaxed dislocation systems is when wc assume that 
the dipole momenta of the dislocation dipoles varies according to a given distribution. 
Let -Pdip(yo) denote the probability distribution of the dipole height yo- Like in the 
monodispcrsc systems studied above (where Pdip(yo) is simply a Dirac delta function), 
the correlation function can be explicitly given by 

dt'ir) = —PdUy)5{x - y - £+(re,t, y)) - ^- (43) 

Pdis -iV 

For the case of the four different types of dipoles the correlation function could be given 
accordingly. Prom (41) and (42) one concludes that 

a'(re.t) ^ -1-53^, (44) 
/3'(rext) «^ -3-16^^ J PaMy'dy. (45) 

3.4. Model of real relaxed dislocation configurations 

As it was mentioned above, relaxed 2D configurations are made up of two basic 
structures, dipoles and walls. In a simple picture, a certain fraction of the dislocations 
forms dipoles with different momenta. Let us denote this ratio by K. The remaining 
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dislocations are distributed in the walls. This means that the correlation function given 
by (43) of the disperse case has to be simply multiplied by K: 

dt~{r) = —PdM6{x -y- 5+(rext, y)) - ^ (46) 
and so for systems containing all four ideal dipole types in equal number 

a'(rext) ^ -1-53^, (47) 

/3'(rext) « -3-16^ J PaMv'dy. (48) 

R 

Csikor and Groma determined the value of K by investigating relaxed configurations 
[23]. They looked at nearest neighbour dislocations, and measured how frequently their 
signs are opposite. According to their result 

= 0.76 ± 0.01 . (49) 

A -Pdip(l/) dipole height distribution was also established by measuring the distance of 
the nearest dislocation of opposite sign for each dislocations. For relaxed configurations 

„ / N 1 f \y\\ . n 0.35 

P^ipiy)^ —ewi-—] with yo~^= 50 



/Pdi; 

was found. With these, from (47) and (48) 

«'(rext)~-1.16-^, (51) 



and 



- (52) 



4. The evolving distribution function 

In the considerations given above the asymptotes of the Fourier transform of the 
distribution function are given. For its real part the result obtained by Csikor and 
Groma given by (12) [23] holds for the loaded case too, while for the complex part (33) 
was derived. By combining (12) and (33) we arrive at 



1 + Cp,,,q' In (M^ + i (G6)^^^pdisg^ In ) , if \q\ - 0, 

exp I -— pdis|g| 1 exp \ \Gh — - — pdis? 1 , if m oo. 



(53) 



The a' and (3' functions, appearing in the above expressions, were calculated in section 
3 for different model dislocation systems. (The function does not have any effect on 
the asymptotes of the resulting distribution function, so it was not given explicitly.) 

As the last step, in this section the asymptotic decay of Pstr('?') and the behaviour 
of its centre part is determined from its Fourier transform given by (53). 
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4.1. Behaviour of the third order restricted moment and the asymptotic decay of the 
distribution function 

Let us divide the |g| — > case expression in (53) into two terms: 

PsHq) = PLAq) + PLA^) (54) 

with 

nli(?):=l + Cpdisg'lnfM') (55) 

and 

PsUl) i (Gbf^^P^isq' In (^^) . (56) 

By analysing (55) Groma and Bako [22] have found that 

1 

-Pstr,i(T) = Cpdisr-pr , if |t| oo. (57) 

r\ 

In order to find Pstr,2, first, we recall the well known fact that for every distribution 
function p differentiable k times the identity 

mfe = i'=(p^)(^)(0) (58) 

holds, where is the A;th order moment of p: 

oo 

™.:=/.VWd.. (59) 

— OO 

and (•)('^) denotes the kth derivative. 

We have to note, that the above relation cannot be applied for P^r,2 since its 3rd 
derivative is infinite at q — 0. It can be seen, however, if we define the function 

P!^2i<l) ■= PLM - Ys^sUm = i (G6)3^^pdisg'ln2, (60) 

unlike P^^ 2) it has a finite 3rd derivative at zero. So, the relation (58) can be applied 
with k = 3 leading to 

T 

lim /(T')'5r2(^')dT' = 3(G'6)V(re.t)Pdisln2. (61) 



T— >0O 

— T 



On the other hand 

T 



/ (r')^:^(r')dr' = J {tT PsuAr'W - ^ | {r'fPsuAr'mdr' 

—T — T — r 

^V3(r)-v^{T/2), (62) 



where the function 



r 

vsir) J (r')^Pstr,2(r')dr', (63) 
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called the third order restricted moment, is introduced. After combining (61) and (62) 
one gets: 

Urn [vsir) - vsir/2)] = 3(G'6)=^a'(Text)Pdis In 2, (64) 

r— >oo 

which is a function equation for V3{t). Its general solution is 

lim vs{t) = 3(G6)3a'(Text)pdis In ( -] , (65) 
where Tq is an arbitrary constant. It follows that 

p . ^ 3(G'&)^a^(rext)Pdis .r I I 

^str,2(r) = —-^ , if \t\ OO. (66) 

Finally, from (54), (57) and (66) one concludes that 

Pstr(r) = Pstr,i(r) + Pstr,2(r) = ^ + ^, if |r| ^ OO, (67) 



where 



U = ^{Gb)'a'{T,^,). (68) 



4-2. The central part of the distribution function 

According to (53), in the limit of |g| oo, due to the external stress, the Fourier 
transform is multiplied by the phase factor exp (i Gb^-^^Y^ PdisQ^ ■ From the identity 

^ y /(a; + Xo)e-'"'''dx = /^(g)e'«^° (69) 

R 

valid for an arbitrary function /, one can speculate that at |r| — > the distribution 
function is shifted with At := Gb^-^^^^ Pdis- However, in (69) the factor e'^^° is required 
for all q, not only in the limit of |g| — > oo. This means, that the value of At given above 
must be considered only an approximation (for details see section 5). To sum up 

Pstr(r) = Pstr,rext=o(r + At), if |t| ^ 0, (70) 

where Pstr,Text=o denotes the stress distribution function of the unstressed case, and 

At ^ Gb^^p^^. (71) 
5. Numerical results 

In this section, we present numerically obtained stress distribution functions and 
compare them with the theoretical predictions discussed above. First, the monodisperse 
case, for which analytical solution was given in section 3, is considered. Afterwards, real 
relaxed 2D configurations are investigated. 
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5.1. The stress distribution function in monodisperse dipole systems 

In section 3.2 the monodisperse dislocation dipole configuration consisting of randomly 
distributed dipoles having the same slip line distance was introduced and analysed. 
Each dipole was deformed separately due to the externally applied stress (see figure 8). 
For symmetry reasons the numbers of the dipoles with the four possible dipole angles 
45°, 135°, 225° and 315° were taken to be equal. 

The reason for analysing the monodisperse dipole system, is the fact that in real 
relaxed dislocation configurations is mainly influenced by the narrow dipoles (see 
section 3.1). In this model system the ci^~ function can be analytically given [see (38), or 
(43) if the dipole height has a certain distribution], which permits to obtain an analytical 
result for the stress distribution function [(41), (42), (67), (68), (70), and (71)]. In this 
section this analytical prediction is validated, proving the correctness of the presented 
deduction. 

A monodisperse dipole configuration can be easily generated numerically according 
to its formal definition in section 3.2. An example is seen in figure 10. In order to 
determine the distribution function of internal stresses, the stress values at the grid 
points of a smooth squared mesh were determined using (1). For the stress field (2) 
was used. (It must be noted that it was numerically verified that if one uses Tind 
corresponding to periodic boundary conditions, it does not affect the results presented 
in this section.) To achieve better statistics, the average of the distribution functions 
obtained on many different monodisperse configurations consisting of 512 positive and 
512 negative dislocations was taken. A typical result can be seen in figure 11. 



At first glance, there is only a little difference between the distribution functions 
in the stressed and unstressed cases. To prove the existence of the additional like 
term in the tail of the distribution, the third order restricted moment, defined by 
(63), was calculated. According to (65), plotting v^^t) as a function of ln(r) must result 




Figure 10. A monodisperse dipole configuration of 1024 dislocations under applied 
external shear stress. The dipole configurations correspond to the equilibrium 
condition (37). 
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a. 
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T/(Gfo) 

Figure 11. The distribution of internal stresses in a monodisperse dislocation dipole 
system subjected to external shear stress Text = (5/-\/1024) • G&^Pdis ~ 0.156 • Gh^p^^s- 
The result was obtained by averaging over 2000 different configurations. Previous result 
for the tail at zero stress level is also plotted [23] . 



in a linear curve. Its slope should be proportional to the external stress. In figure 12 
this method is applied for different external stress values. The existence of the term 
is clearly confirmed. The linear dependence of the coefficient a' (Text) on the external 
stress Text is also verified in figure 13. By fitting a straight line one obtains 

U = -2.36iGb)\,t. (72) 
From (68) we get 

a'(rext) = -1-57^ (73) 

meaning there is only about a 2.5% difference between the numerically and theoretically 
obtained coefficient of the r~'^ tail [see (41)]. 

It was derived theoretically that due to the external stress Text the central 
Lorentzian-like part of the stress distribution function is shifted with At, which is 
proportional to Text [see (70) and (71)]. The shifting can be seen in figure 14. The At 
values corresponding to different external stresses are plotted in figure 15. The linear 
stress dependence of the shifting is clearly fulfilled. Fitting a straight line yields 

/5'(rext) = -2.67^yo' (74) 

for Ho = 0.32 ■ Pdis'^. There is about a 20% difference between this numerical result and 
the theoretical prediction given by (42) and (71). As we mentioned already in section 
4.2, (71) serves only as an estimate for the shifting, therefore, the error of 20% means 
that numerical results are in a quite good agreement with the theoretical prediction. 

5.2. The stress distribution function in relaxed systems 



In this section we repeat the previous calculations for real relaxed 2D configurations 
obtained by discrete dislocation dynamics simulations. Initially 64 positive and 64 
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1.5-10^ 




Figure 12. The vs restricted moments calculated from the distribution functions 
according to (63). By plotting W3(t) as a function of ln(r) the resulting linear curve 
proves the t~'* like term in the tail of the distribution function [see (65)]. The 
method was repeated for different applied external stresses Text (measured in Gh^pdis 
dimensionless units). The results were obtained by averaging over 2000 different 
configurations. 



Numerical data 




-0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 

Figure 13. The values of U obtained by fitting to v-^It) curves in figure 12 for 
different external stresses Text [see (65) and (68)]. The linear dependence obtained is 
in agreement with the theory [see (41)]. 

negative dislocations were placed randomly in a square domain, then they were let to 
relax under applied external shear stress with the conventional overdamped dynamics 
[1, 7, 38]. To emulate infinite medium, periodic boundary conditions were used (for 
details about the stress field generated by a dislocation under these conditions see [39]). 
In order to reduce computational time, annihilation was introduced for dislocations of 
opposite signs if their distance became less then 0.02pJ°'^. This affected in average only 
the 5% of starting dislocations. The stress distribution function was determined in the 
same way described in the previous section. Again, it was found that the distribution 
function is not infiuenced by the boundary conditions imposed for the stress. To 
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Figure 14. The central part of the stress distribution function of a nionodisperse 
dislocation dipole system under different applied external shear stresses. The result was 
obtained by averaging over 2000 configurations. (The external stress Text is measured 
in Gb^pdis dimensionless units.) 
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-0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 

Figure 15. The shift of the Lorcntzian describing the central part of the stress 
distribution function of a monodisperse dislocation dipole system under different 
applied external shear stresses. Again, the distribution function is the average of 
2000 different realizations. 

obtain acceptable results 2000 parallel simulations were carried out and averaging was 
performed over the stress distribution functions. 

Like in the previous section, in order to prove the asymptote (67), in figure 16 the 
third order restricted moment vs^t) was plotted against ln(r) [see (65)]. Due to the 
huge computational demand of the calculations, simulations were only performed with 
two different external stresses. At Text = for symmetry reasons /'str(^) — -Pstr( — 
Therefore v^lr) must vanish. For opposite, non-zero external stresses, again for 
symmetry reasons, -Pstr,rcxt (''") = -^str,- Text ("''")• This implies that 

^3,rext(^) = -^3,-rext(^)- (75) 





24 



The values of U were determined according to (65) and (68). The results are plotted in 
figure 17 together with the theoretical prediction given by (51). It can be seen that the 
coefficient U indeed depends linearly on the external stress. By fitting a straight line to 
the data points we get 

U = -0.85(G6)Vext. (76) 
According to (68) 

«'(rext) = -0.57^. (77) 
Go 

This means that the measured slope of the ?7(rcxt) function is about half of the one 
calculated from the dipole approximation. 
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Figure 16. The ?;3 restricted moments calculated from the distribution functions 
according to (63). By plotting W3(t) as a function of ln(T) the resulting linear curve 
proves the like term in the tail of the distribution function [see (65)]. The method 
was repeated for two different applied external stresses Text (measured in Gfe^pdis 
dimensionless units). The results were obtained by averaging over 2000 different 
relaxed configurations. 

The difference observed can be attributed to the fact that in the considerations 
explained above it was assumed that the relaxed system is mainly built up from small 
dislocation dipoles deformed only due to the external stress, i.e. the dipole-dipole 
interaction was neglected. In real systems, however, this is not the case. If one looks 
at figure 7(b), beside dipoles, dislocation multipoles can also be observed in a large 
number. As it is indicated in figure 18, a dislocation multipole is always 'harder' than a 
dipole, which formally means e^(rext,yo) < ^^(^ext, l/o) (see figure 18 for the definition 
of e^). For a small external stress, however, e^(rext,2/o) oc Text still holds. So, the 
asymptotic decay of the distribution function remains the same, only the coefficient in 
(51) is smaller. 

Concerning the shift of the central part of the distribution function (see figure 19), 
it can be seen in figure 20 that the shift Ar is proportional to the external stress with 

/?'(rext) = (78) 

Pdis Gh 
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Figure 17. The U values obtained by fitting to V3{t) curves in figure 16 for different 
external stresses [see (65)]. The obtained linear dependence is in agreement with the 
theory [see (51)], however its coefficient is different (the explanation is discussed in 
the text). For Text < values the Ur„^^ = —U-r^^^ relation was used (see the text for 
details). 




yo 



Figure 18. The deformation of a certain dislocation multipole under Text external 
shear stress. The change of the x coordinate of the negative dislocation is denoted 
by e+(roxt,yo) which is here half of the similar deformation of a single dipole: 
ei(T-cxt,yo) = e+(Toxt,|/o)/2. 



6. Conclusions 

A detailed analysis of the distribution of internal shear stresses in a 2D dislocation 
system subjected to external shear stress was presented. The results can be summarised 
as follows: 

(i) It was shown theoretically that due to applied stress the initially symmetric stress 
distribution function becomes asymmetric: 

(a) A term proportional to l/(r|r|^) is added to the 1/|t|'^ like tail of the 
distribution function. The coefficient of the extra term is proportional to the 
stress applied. 
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Figure 19. The central part of the stress distribution function of relaxed dislocation 
caused by external shear stresses. The result was obtained by averaging over 2000 
different configurations. 
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Figure 20. The shift of the Lorentzian describing the central part of the stress 
distribution function of a relaxed dislocation dipole system under different applied 
external shear stresses. 

(b) The central Lorentzian like part of the stress distribution function [23] is shifted 
with a value proportional to the applied stress. 

(ii) For the monodisperse dislocation dipole system the theoretical predictions were 
proved by the numerical calculation of the stress distribution functions. 

(iii) The distribution function was determined for 2D relaxed dislocation configurations 
generated by discrete dislocation dynamics simulations at different applied external 
stresses, too. It was found that one must step beyond the dipole approximation 
because of presence of dislocation multipoles. However, they only make the material 
'harder', in the sense that larger external stress is needed to achieve the same change 
in the distribution function as in monodisperse dipole systems. 
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